A statistical study of multiply-imaged systems in the lensing cluster Abell 68.^'^ 

Johan Richard^, Jean-Paul Kneib'^'^, Eric JuUo^'^, Giovanni Covone^, Marceau Limousin^, Richard EUis^, 
Daniel Stark^, Kevin Bundy^, Oliver Czoske^, Harald Ebeling^ and Genevieve Soucail^" 

ABSTRACT 

We have carried out an extensive spectroscopic survey with the Keck and VLT telescopes, 
targeting lensed galaxies in the background of the massive cluster Abell 68. Spectroscopic mea- 
surements are obtained for 26 lensed images, including a distant galaxy at z = 5.4. Redshifts 
have been determined for 5 out of 7 multiply-image systems. Through a careful modeling of 
the mass distribution in the strongly-lensed regime, we derive a mass estimate of S.SxlO^^M© 
within 500 kpc. Our mass model is then used to constrain the redshift distribution of the re- 
maining multiply- imaged and singly-imaged sources. This enables us to examine the physical 
properties for a subsample of 7 Lyman-a emitters at 1.7 < z < 5.5, whose unlensed luminosities 
of ~ 10"*^ ergss^^ are fainter than similar objects found in blank fields. Of particular interest is 
an extended Lyman-a emission region surrounding a highly magnified source at z = 2.6, detected 
in VIMOS Integral Field Spectroscopy data. The physical scale of the most distant lensed source 
at z = 5.4 is very small (< 300 pc), similar to the lensed z ~ 5.6 emitter reported by Ellis 
et al. (2001) in Abell 2218. New photometric data available for Abell 2218 allow for a direct 
comparison between these two unique objects. Our survey illustrates the practicality of using 
lensing clusters to probe the faint end of the 2^2 — 5 Lyman-a luminosity function in a manner 
that is complementary to blank field narrow-band surveys. 

Subject headings: cosmology: observations — galaxies: clusters: individual (A68) — gravitational 
lensing — galaxies: high redshift 
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Introduction 



The central regions of massive galaxy clusters act as powerful gravitational telescopes, magnifying the 
light from background galaxies via the effect of strong lensing. Such magnifications can attain typical values 
of 1 to 3 magnitudes in concentrated cluster cores, enabling the detection of intrinsically fainter sources than 
in unlensed field surveys. The detailed study of low luminosity galaxies at z > 2, where the major fraction 
of star-formation activity is thought to occur, is an intere sting but poorly-understood topic. S uch galaxies 



can either be found through their Lyman-a emission (e.g.. iFranx et al 



their ultraviolet continuum fluxes via the Lyman break techniques (|Kneib et al.ll2004 ; 



1997 



Santos et al. 2004). or through 



Richard et al. 200' 



A prerequisite for strong lensing studies of intrinsic ally faint galaxies at high redshift is an accurate 
meas urement of the projected mass distribution in the lens (jKneib et al iGavazzi et al l200,4 I Sand et al 



20051) . 



ed sources of 



known redshift. 



2006 



Broadhurst et al. 


2005; 


Halkola et al. 


ems, or Abeh 2218 (Ebbels et al. 1996; 



Kneib et al.lll996| ) have sufficient constraints to permit precise modelling of each individual dark matter 
clump. 

Spectroscopic searches for Lyman-a emitters (LAEs) at high redshift usually have a better line flux 
sensitivity and span a larger redshift range (Az 4) than those of wide-field narrow-band surveys. This 
gain in sensitivity is even larger in strong lensing applications. Lensed spectroscopic surveys may also be 
sensitive to sources with emission lines with an equivalent wi dth W < 20 A , smaller than those in narrow- 
band surveys (e.g., Fvnbo et al. 2003t Shimasaku et al. 2006). An additional complication in narrow-band 
surveys is how interlopers are treated; confirmatory spectroscopy is usually necessary. By contrast, in lensed 
surveys, the geometrical configuration of multiply- i maged systems c an reliably distinguish between high 
redshift objects and low redshift interlopers (see e.g., lEllis et ai.J|20Qll). 

As our surveys expand, a variety of types of emission line galaxies are being discovered. Of particular 
interest are the extended Lyman-a emission sources which have been mainly discovered in re gions of signifi- 



cant o verdensity through deep narrow-band imaging ISteidel et allboool : IPrancis et allboOll ). 



Matsuda et al 



(j2004[ ) have identified a large number of such giant Lyman-a blobs (with a typical size > 50 kpc) in a 34' x 27' 
field of view, demonstrating the existence of a continuous distribution. The origin of the extended Lyman-a 
emission in such radio-quiet sources may be explained by gas inflow during the early stages of galaxy forma- 
tion: large amounts of hydrogen collapsing into the dark matter potential well will cool through Lyman-a 
radiation. Giant Lyman-a blobs may thus be the progenitors of very massive galaxies in the local Universe. 
A key issue is whether the same process is seen to occur in lower-mass objects. A route to addressing 
this question is to examine the nature of smaller extended Lyman-a sources, either by long-slit or Integral 
Field Spectroscopy (IFS). This identification is more easily accomplished in strongly- lensed sources where 
magnification will stretch the observed physical scales. 

The spatial magnification associated with lensing can als o be used to yield physical sizes for the most 
distant sources. Using strong lensing in the cluster Abell 2218. lEllis et al.l (j200ll ) located a remarkably small 
source at 2;=5.6 where the combination of the Lya emission line flux density and the weak stellar continuum 
were used to deduce a young age and modest stellar mass (~ 10^~^ -^o) consistent, perhaps, with a forming 
globular cluster. Further surveys are required to evaluate whether such systems are common at z ~6. 

The major drawback arising from the study of lensed sources located through studies of individual 
clusters is, of course, the significant cosmic variance that is associated with the small volumes being probed. 
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Compared to field surveys, any statistical inferences on the abundances of various classes of populations 
may be much more uncertain, even granting fainter sources are probed. To overcome this limitation, an 
effective survey would have to be conducted through a large sample (~20-40) of lensing clusters, each wit h 
reliable mass models based on the spectroscopic study of many multiply-imaged systems ( Kneib et al. 1996[ ). 



Fortunately, the construction of such a sample of well-mapped clusters is now a realistic proposition. Several 
Hubble Space Telescope (HST) snapshot imaging surveys of X-ray luminous clusters are now underway with 
associated ground-based spectroscopy, such as the MAssive Clusters Survey (MACS, GO#10491, P.I.: H. 
Ebeling) and the Local Cluster Substructure Survey (LoCuSS, GO#10881, P.I.: G. Smith). 

The purpose of this paper is to illustrate the promise of such surveys by examining spectroscopically the 
rich population of lensed sources located in the lensing cluster Abell 68 (a=00:37:06.81 15=-|-09:09:24.0 J2000, 
z = 0.255), one of the most X-ray luminous clu sters (Lx ~ 8.4 ±2.3 x lO^'^'ergs"^, 0.1-2.4 keV) in the X-ray 



Brightest Abell- type Clusters samp le (XBACS jEbeling et al.l (jl996[ )). Strong lensing in this cluster has been 



previously studied by Smith et al. ( 2005h . hereafter S05, as part of a survey of 10 X-ray luminous galaxy 



clusters at z ~ 0.2. Smith et al identified a list of potential multiple-image systems, a few of which were 
confirmed spectroscopically. Here we significantly extend this work by securing the redshifts of new multiple- 
image systems, many of which are strongly-lensed Lyman-a emitters at z > 2. The combination of a large 
magnification factor, high-resolution HST imaging and broad-band photometry enables us to demonstrate 
the value of studying the physical properties of these faint emitters, such as their star formation rates, 
intrinsic scales and stellar masses. The paper is intended to illustrate the significant promise of continuing 
such spectroscopic work with the larger samples of clusters now being surveyed with HST. 

The paper is organized as follows. In Section 2, we describe the various observations and the reduction 
of the spectroscopic data. We present in Section 3 the strong- lensing constraints, in the light of the redshifts 
and identification of new multiply-imaged systems. Section 4 presents a mass model of the cluster from 
which the source magnifications are deduced. The physical properties of the various categories of high 
redshift LAEs are presented in Section 5 and the implications are discussed in the context of the limitations 
of blank field surveys in Section 6. We summarize our conclusions in Section 7. 

Throughout this paper, we adopt the following cosmology: a flat A-dominated Universe with the values 
r^A = 0.7, = 0.3 = 0.04 5 and Hq = 70 kms"^ Mpc~^ All magnitudes given in the paper are quoted 



in the AB system (jOkd Il974l ). The correction values Cab between AB and Vega photometric systems, 
defined as mAB = wvega + Cab, are reported in Table [1] for each filter. At the redshift z — 0.255 of the 
cluster, the angular diameter distance is 3.9 kpc arcsec"^. 



2. Observations and Data Reduction 

We present in this section the photometric and spectroscopic datasets used to assemble our catalog. 
High resolution images are crucial for the morphological identification of multiple image systems and the 
precise astrometric position of the sources studied here, whereas multicolor images are used to estimate their 
spectral energy distributions. Redshift and emission lines measurements for individual objects were obtained 
during subsequent spectroscopic observations. These included multi-object spectroscopy of multiply-imaged 
candidates, as well as systematic long-slit searches in the central regions of the cluster. Figure [1] shows the 
location of the main spectroscopic settings in the cluster field. 
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2.1. Imaging data 

A considerable body of multi-wavelengtii data exists in the field around Abell 68, including high reso- 
lution HST imaging. The main characteristics of the dataset used in this study are summarized in Table [1] 
3 X 2.5 ks of integration time with the Wide Field Planetary Camera (WFPC2) was obtained during Cycle 
8 in the i?-F702W band, as part of HST program #8249 (PI : J. P. Kneib). Observations were carried out 
in low sky mode, and a 1.0" dithering pattern was used between each exposure. Details on the reduction of 
these data are given in SOS. 



Recognition of faint multiply- imaged systems in the vicinity of the cluster core is hindered by the dom- 
inant stellar halo of the Brightest Cluster Galaxy (BCG). To overcome this, we fitted and subtracted from 
the HST image a model representation of the surface brightness distribution using the IRAF task ellipse. 
Both the position angle and ellipticity were allowed to vary as a function of the semimajor axis in the fitted 
elliptical isophotes, as well as the isophote centroid in the central part. This procedure was found to give 
satisfactory residuals at the center (Figure [3]). 



Associated optical images in B, R, I have been obtained on UT 1999 November 19 using the CFH12k 
camera at CFHT. These sample a field of 42'x28' at a 0.205" pixel scale. The total exposure times are 8.1, 
7.2 and 3.6 k s in the B, R and J band, respectively. T he data was reduced using procedures similar to those 



described bv lCzoskel (|2002h and lBardeau et al.l |2005h 



At longer wavelengths, Abell 68 has been observed at the Very Large Telescope using the FOcal Reducer 
/ low dispersion Spectrograph (FORS2/UT4) in the z-band on UT 2002 October 06, and the Infrared 
Spectrometer And Array Camera (ISAAC/UTl) in the J and H bands on UT 2002 September 29 . The 
field of view of the F0RS2 image is 7.2 x 7.2 arcmins after dithering, with a pixel size of 0.252", and we 
used 80 dithered exposures of 120 s. The field of view of the ISAAC images is about 2.5 x 2.5 arcmins after 
dithering, with a pixel size of 0.148", the subintegration x integration times of the dithered exposures were 6 
X 35 s and 10 x 12 s in the J and H bands, respec tively. All these data have been reduced using procedures 
similar to those described by Richard et al. ( 2006h . 



Instrument 


Filter 


Exposure Time 


Pixel Size 


Depth 


Cab 


Seeing 






(ks) 




AB mag. 


mag. 


II 


CFH12k 


B 


8.1 


0.206 


27.4 


-0.066 


1.11 


CFH12k 


R 


7.2 


0.206 


27.2 


0.246 


0.67 


WFPC2 


R702W 


7.5 


0.1 


28.0 


0.299 


0.17 


CFH12k 


I 


3.6 


0.206 


26.5 


0.462 


0.58 


F0RS2 


z 


9.6 


0.252 


26.5 


0.554 


0.71 


ISAAC 


J 


6.48 


0.148 


26.2 


0.945 


0.48 


ISAAC 


H 


7.12 


0.148 


26.3 


1.412 


0.48 



Table 1: Properties of the photometric dataset: from left to right: instrument and filter names, total integra- 
tion time, pixel size, photometric depth (defined as 4 pixels above 3 cr, where a stands for the typical local 
background noise), photometric correction Cab between AB and Vega systems, seeing measured on bright 
unsaturated stars. 
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Fig. 1. — Composite CFH12k-i3 i?J color image of the field of view around the eenter of Abell 68. We 
overplot the redshift measurements obtained for galaxies located in the background of the cluster (black 
labels). Red circles represent cluster members confirmed with spectroscopy. We delineate the imprints of 
the HST/WFPC2 (black polygon) and the VIMOS/IFU (blue square) fields, as well as the spatial coverage 
of the different LRIS long-slit configurations (red rectangles). 
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2.2. Keck Multi-Object Spectroscopy 



The Low Resolution Imaging Spectrog raph (LRIS;[0keei3 Hiii)) on Keck I has been used in multi 



object (MOS) mode during two observing runs, in order to target background galaxies and multiply-imaged 
candidates selected on the basis of morphology and colors. 

On UT 2001 August 4, 4 exposures of 1.8 ks were acquired with a 31 slits-mask, using a 300 1 mm~^ 
grism blazed at 5000 A in the single red channel of the camera, which covers the approximate range 5500- 
9900 A at the dispersion of 2.5 A per pixel. The average seeing was ~ 1.0". The night was photometric and 
spectrophotometric standard stars were used for flux calibration. 

On UT 2002 November 30, a 32 slits-mask was used during 3 x 2.4 ks of integration time. A 6800 A 
dichroic separated the red channel of the instrument, equipped with a 600 1 mm^^ grating blazed at 7500 
A, from the blue channel equipped with a 400 1 mm~^ blazed at 3400 A. The whole setting covers the 
wavelength range 3500-9500 A with dispersions of 1.28 and 1.09 A per pixel in the red and blue channels, 
respectively. Despite good seeing conditions 0.8"), the night was not photometric and no standard stars 
were observed. 

These datasets were reduced using standard IRAF procedures for bias removal, flat-fielding, wavelength 
and flux calibration, sky subtraction and extraction of the one-dimensional spectra. 



2.3. Keck long-slit spectroscopy 

2.3.1. Optical spectroscopy 
Abell 68 was observed on UT 2002 Sept ember 11 with LRI S, in the course of a survey targeting low- 



luminosity Lyman-a sources at high redshift (jSantos et al.ll2004^ . A 175"-long and l"-wide slit was used to 



map the high magnification regions of a sample of z 0.2 lensing clusters. In the case of Abell 68, 6 adjacent 
slit settings scanned the theoretical location of the critical li nes at z ~ 5, with 2x1000 s of integration time 



at each position. The reduction of these data is detailed in (jSantos et al.ll20041) 



In addition to the detection of high redshift sources through their Lyman-a emission, this blind spec- 
troscopic survey provided secure redshifts for a number of lensed background galaxies serendipitously falling 
into the long slit. 

On UT 2003 August 26, a single long-slit LRIS position was aligned on two components of a triply- 
imaged system, discovered as R-dropouts by a RIz color-color selection technique (see Sect. [3|). A 300 1 
mm~^ grism blazed at 5000 A and a 600 1 mm~^ grating blazed at l/.tm were used in the blue and red 
channels of the instrument, both lightpaths being separated by a dichroic at 6800 A. Two exposures of 1.2 
ks were acquired at this position, with a 5" dithering offset along the slit. 

Finally, an additional LRIS long-slit integration of 3 x 1.2 ks was acquired on UT 2005 November 29 
with a 600 1 mm~^ grism blazed at 4000 A, a 400 1 mm ^ grating blazed at 8500 A, a 5600 A dichroic, and 
5" dithering offsets. 
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2.3.2. Near Infrared Spectroscopy 



Abell 68 was ob served on UT 2005 October 13 using the Near InfraRed SPECtrograph (NIRSPEC, 



McLean et al.l (|l998r )) on Keck II, during a spectroscopic survey of t he critical lines of lensing clusters 
similar to the LRIS survey described before, but at longer wavelengths |stark fc Elhslbood) . A 42"x 0.76" 
long slit was used at 2 adjacent slit positions, with 9 x 600 s of integration time on each of them. The 
spectra were r educed using IDL routines, following optimal spectroscopic reductions techniques presented in 
KelsonI (|2003h . More details are presented in a forthcoming paper (Stark et al. 2006b, in press). 



2.4. VLT - Integral Field Spectroscopy 



Abell 68 was ob served on UT 2004 August 12 using the Visible Multi-Object Spectrograph (VIMOS, 
Le Fevre et al.l (j2003t ll on VLT/UT3 in low-resolution (LR-Blue grism) Intregral Field Spectroscopy mode, as 
part of a survey targeting the central regions of an intermediate-redshift galaxy cluster sample (073.A-0774, 
PI: Soucail). The 54" x 54" field of view of the Integral Field Unit (IFU, composed of 6400 fibers, splitted 
in 4 quadrants feeding the 4 VIMOS CCDs) was centered on the cD galaxy. In the given configuration, 
the spectral resolution is about 200 and the diameter of the fibers is 0.66", covering the wavelength range 
3900-6800 A with a dispersion of 5.355 A per pixel. 2 x 2.4 ks of integration time have been acquired without 
dithering. 

The se 3D spectroscopic da ta have been reduced using the Vimos Interactive Pipeline Graphical Interface 



(VIPGI. IScodeggio et al.ll2005r l^l. Before building the data cube for each exposure, every step in the data 
reduction was perforin ed on a single quadrant basis. After bias subtraction and cosmic ray removal (see 



Zanichelli et~ai1 |2005h for a description of the algorithm), the spectra were traced on the CCDs with the 
help of the high S/N spectra of a continuum lamp. Following wavelength calibration, inhomogeneities in 
fiber efficiencies were corrected by measuring the counts in the 5577 A sky emission line after subtracting the 
contribution from a galaxy spectrum where fibers cover a ga laxy position. The flu x-calibration is applied by 
using observations of a standard star in each quadrant. See Covone et al. ( 20061 ). who present similar data 
on Abell 2667, for a detailed description of the procedures. 



2.5. Redshift measurements 



We attempted to measure the redshift of all individual objects falling in the slits that revealed a dis- 
cernable continuum or possible emission lines. To obtain an accurate redshift measurement for foreground 
galaxies, cluster members and other bright objects, we applied the IRAF task xcsao from the Radial Velocity 
package RVS AO (jKurtz et al. 1992 ) on all extracted spectra. This pro cedure uses a cross-correlation method 
based on spectral templates ( Coleman et al. 198d Kinnev et al. 1996h to estimate the redshift and the corre- 
sponding redshift error. For the remaining objects in the spectroscopic catalog, the redshift measurement is 
based on the wavelength at the peak of the brightest emission line detected. In the latter case, we estimated 
the redshift error from the spectral dispersion. Additional uncertainties generated by the accuracy of the 
relative and absolute wavelength calibrations, of about 0.8 and 1.5 A respectively for the LRIS data, were 



^ VIPGI has been developed within the VIRMOS Consortium. 
[http:/ /cosmos. mi. iasf.cnr.it/pandora/vipgi. html. 



For more information, see 
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ID 


RA 


Dec 


R 


z 


C 


Features 


M 


Wo 




00: 


09: 


[mag] 








[mag] 


[A] 


















(Lya) 


1 (C15a) 


37:04.297 


09:43.40 


26.00 ±0.18 


5.421 ± 0.0021 


9 


Lya 


2.74±0.08 


53 ± 16*^ 


2 (C15b) 


37:04.861 


09:51.78 


26.52 ±0.22 


5.421 ± 0.0021 


9 


Lya 


2.89±0.07 




3 (C26) 


37:08.960 


09:06.90 


27.24 ±0.25 


3.677 ±0.0022 


9 


Lya 


2.06±0.03 


107 ±8.2 


4 (C23a) 


37:07.902 


09:29.80 


26.44 ±0.18 


3.135 ± 0.0021 


9 


Lya 


2.47±0.07 


25.9 ±5.2 


5 (C25) 


37:06.506 


10:16.70 


24.40 ± 0.07 


2.6930 ± 0.0021 


9 


Lya 


1.12±0.07 


10.4 ±5.1 


6 (C20c) 


37:05.405 


09:59.14 


25.15 ±0.09 


2.6890 ± 0.0020 


9 


Lya 


3.61±0.09 


30.4 ±6.3 


7 (C4) 


37:07.657 


09:05.90 


23.31 ±0.04 


2.6280 ± 0.0021 


9 


Lya 


4.15±0.16 


42.4 ±2.3 


8 (C27) 


37:04.906 


10:30.20 


22.81 ±0.03 


1.7546 ± 0.0021 


4 


Lya 


1.72±0.10 


2.4 ±0.5 


















([On]) 


9 (Cla) 


37:06.207 


09:17.49 


24.02 ± 0.04 


1.5836 ±0.0011 


4 


[Oiii], H/3 


2.52±0.06 


< 50 


10 (C12) 


37:04.930 


10:21.40 


21.67 ±0.02 


1.0171 ±0.0011 


3 


[On], Ha 


1.63±0.06 


10.5 ±1.4 


11 (C7) 


37:05.073 


10:04.90 


22.78 ± 0.03 


0.8610 ± 0.0007 


3 


[On] 


1.59±0.02 


21.3 ±5.9 


12 (C8) 


37:03.700 


09:54.10 


23.08 ± 0.02 


0.8609 ± 0.0008 


4 


[On] 


1.37±0.03 


112.8 ±5.1 


13 (C24) 


37:05.900 


09:59.70 


23.64 ± 0.02 


0.8152 ±0.0007 


1 


K, H, H7 


1.33±0.01 


N/A 


14 


37:04.352 


07:39.60 


21.91 ± 0.03 


0.6386 ±0.0007 


3 


[On] 


0.10 


14.1 ±4.5 


15 


37:08.541 


08:01.30 


22.49 ±0.03 


0.6281 ±0.0006 


3 


K, H 


0.20 


< 10 


16 


37:07.100 


08:23.10 


21.44 ± 0.02 


0.6275 ± 0.0006 


3 


K, H 


0.28 


14.6 ±2.2 


17 (C14) 


37:08.534 


09:14.10 


22.02 ±0.02 


0.6234 ± 0.0007 


2 


K, H 


0.93±0.01 


16.5 ± 7.4 


18 


37:02.707 


08:19.18 


20.52 ± 0.02 


0.6225 ± 0.0006 


3 


K, H 


0.20 


6.7 ±1.3 


19 


37:05.505 


09:24.24 


22.71 ±0.02 


0.6224 ± 0.0006 


4 


[On], [Oni] 


1.83 


25.3 ±6.5 


20 


37:01.890 


07:27.70 


21.40 ± 0.03 


0.6195 ± 0.0005 


4 


Mgii, [Oiii] 


0.10 


N/A 


21 


37:05.000 


07:50.60 


20.33 ± 0.02 


0.6177 ±0.0007 


4 


[On] 


0.10 


10.7 ±1.7 


22 


37:01.400 


05:55.20 


22.10 ± 0.02 


0.5944 ± 0.0005 


4 


H/3, [Oiij] 


0.0 


N/A 


23 


37:08.620 


08:51.40 


23.29 ± 0.04 


0.5817 ±0.0006 


4 


H/3, [Oiii] 


0.51 


< 10 


24 


37:06.410 


09:50.65 


24.66 ± 0.05 


0.4941 ± 0.0006 


4 


H/3, [Oin] 


0.75 


< 10 


25 


37:05.670 


10:04.90 


24.00 ±0.04 


0..3958 ± 0.0007 


4 


Ha, [Oiii], [On] 


0.51 


24.5 ±2.2 


26 


36:57.170 


07:00.80 


21.03 ± 0.01 


0.3693 ± 0.0005 


4 


[On], [Oin] 


0.0 


10.2 ±2.1 



Table 2: Spectroscopic catalog of lensed background galaxies. From left to right: astrometric position (J2000), 
i?-band magnitude in AB system (WFPC/F702W, or CFH12k (in italic)), spectroscopic redshift, redshift 
confidence class (see text for details), main spectroscopic features, magnification factor (in magnitudes), 
rest-frame equivalent width of Lyman-a (top half) or [On] emission lines (bottom half). 



"Measured on the averaged spectrum of CI 5a and CI 5b 
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^700 7750 7800 7B50 
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Lya Sill 01 jCII 







5400 5600 5800 6000 5200 

Wavelength (Ansgtroms) 



Lya Sill CI Ol/Sill Cll 



z=5,421 / 
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1^700 7750 7800 7850 7900 

Wavelength (Ansgtroms) 
Lya Sill 

z = 3.135 




4800 5000 5200 5400 

Wavelength (Ansgtroms) 





4400 4600 4800 5000 

Wavelength (Ansgtroms) 



;Lya iSill iiSill Cll liSilV 
z=1.7546:: i 



3400 3600 3800 4000 4200 4400 

Wavelength (Ansgtroms) 





= 1 .0171 




[oil] 


He 


H, - 


f 
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Fig. 2. — 10" X 10" WFPC/F702W zoomed images with the sht location and extracted spectra for all 
individual z > I sources. Dotted lines outline the emission and absorption features identified in each 
spectrum. Upper red spectra have been smoothed using a cr = 3 pixels gaussian, and shifted in the vertical 
direction for clarity. 
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quadratically added to the previous estimates to yield the final redshift errors. 



A confidence cla ss C, ranging from 1 to 4 was assigned to each individual redshift measurement according 
to the prescription of Le Fevre et al. (1995): this corresponds to a probability level for a correct identification 
of 50, 75, 95 and 100%, respectively. A specific value of 9 is used when only a single secure spectral feature 
is seen in emission. The identification of most z > 2 objects in the catalog is based on only a single emission 
line interpreted as Lyman-a and a confidence class of 9. However, the constraints provided by the lensing 
configuration in case of multiple images (see Sect. [3]) enable us to strengthen most of these interpretations. 

Multicolor photometry was performed for a ll sources identified in the WFPC/F702W band and CFH12k- 
R band images, using the SExtractor software ([Bertin fc Arnoutslll996[ ). Total magnitudes and errors were 
measured on the original images without any resampling or convolution. 

The final spectroscopic catalog of lensed galaxies is presented in Table [51 

Figure [1] displays the location of all sources in the catalog and the spectroscopic configuration of each 
dataset. A red circle marks the location of 44 cluster members from our sample, for which we measured a 
spectroscopic redshift. Long-slit and IFU pointings were mainly focused on the central region of the cluster 
(within '--^SO"), in the vicinity of the critical curves, whereas MOS surveys with LRIS probed a wider field 
of view of 6 X 8 arcmin. 

The redshift distribution of lensed background galaxies appears to be highly correlated. Three sources 
are located in the 2.63 < z < 2.69 range, two sources have z ~ 0.86, and CI system has a. z ^ 1.6 redshift 
similar to CO ISmith et^ |2002b[) . Even more prominent is a group of 8 sources having 0.62 < z < 0.64, 
which lie predominantly at the south of the cluster center, and exhibit similar colors in the composite CFH12k 
image (Figure [1]). 



3. Catalog of Multiply-Imaged Systems 

In this section we update the catalog of multiply-imaged systems in the field of Abell 68, taking into 
account both our new spectroscopic measures and systems without spectroscopy located on the basis of 
their geometrical location and similar colors. For each system, we present a close-up view from the HST-i? 
image in the bottom panels of Figure [3l and summarize the position, photometry, shape parameters and 
magnification in Tabled 



3.1. Morphological Identification 

Multiply- imaged systems can be identified via a visual inspection of morphologies in the central region 
of the HST-i?702i4' image (Fig. [3]) , examined in combination with the broad-band photometric catalogs. 
Our starting point is a detailed study of candidates identified to ^jo2 ^ 25 mag arcsec"^ presented by S05. 
This work revealed 3 main systems (COa,b,c; Cla,b,c; C2a,b) as well as 20 other possible multiply-imaged 
candidates (C3 to C22). We adopt the same nomenclature, extending to include the new systems presented 
here. 

The three images identified as C15, C16 and C17 have been selected as i?-band dropouts on the basis 
of color-color diagrams combining R, I and z-band filters. Indeed, they are very faint in the HST image 
(i?702 ^ 25.7), undetected at shorter wavelengths in the B or R band with CFH12k with a combination 
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Table 3: Properties of the multiple-image systems. In addition to the 6 secured systems, CIO is included as 
a potential candidate (see Sect. 13. ip . From left to right: identification, astrometric position (J2000), ellipse 
shape parameters (a, b, 9), i?702VF magnitude, measured redshift, magnification factor (in magnitudes) 
derived from the mass model. Bracketed values are also predictions inferred from the mass model. 
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Fig. 3. — (Top) Enlarged region of the i?-F702W image showing the location and morphology of the multiple- 
image systems identified in the field of Abell 68. We subtracted the contribution from the brightest cluster 
galaxy (BCG) to assist in the identification of underlying objects, such as the radial arc associated with CO. 
The thin/thick black curves show the location of the critical lines at z = 1.6 and z — 5.42 respectively, as 
inferred from the mass model (§4). We overplot in red new spectroscopic redshifts of background sources. 
(Bottom) Close-up on each component of the multiple systems discussed in the text. 
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of red {R — I)ab ^ 1-7 and blue (/ — z)ab ^ —0.5 colors, as measured with aperture photometry on 
the seeing-matched images. Such a spectral energy distribution (SED) is characteristic of galaxies in the 
redshift range 5 < z < 6, where the Lyman-a break in the spectral continuum occurs between the R and / 
bands. The preliminary mass model (§4) shows that the location of these three images is compatible with a 
multiply-imaged source, thereby strengthening the high redshift interpretation. We designate this system by 
C15 (with three images C15a,b,c that correspond to C15/16/17) and confirmed the components C15a and 
C15b lie at the same redshift using LRIS spectroscopy (see Section [3T2|) . 



The arclet identified as C20 was serendipitously covered during the same LRIS long-slit observation. 
We interpret the extended emission seen in the spectrum (Sect. 13. 2p as Lyman-a at z = 2.689. This is also 
supported by the mass model, which predicts 2 detected counter-images with similar optical colors. We refer 
to this system as C20a,b,c. 

Next to C20a and C20b images, a fainter extended arc C6 was identified by SOS as a possible multiply- 
imaged system at a similar redshift. This arc can be split into 3 components C6a,b,c. 

Finally, we un covered a new len sed image, C23a, close to the cluster center during the critical line 



survey using LRIS (jSantos et al.ll2004( l. Its spectrum is consistent with a high redshift source dominated by 
its Lyman-a emission. Again, the cluster mass model predicts the position of a second faint component for 
this system (C23b), detected on the HST image. 



3.2. Redshift Constraints 

We now present the new spectroscopic redshifts obtained for 6 multiply-imaged candidates, as well as 
constraints implied for the remaining multiply- imaged systems C2 and CIO, from the updated mass model 
(Sect. H. 



CO source 

This source is a strongly-lensed (p 3.2 mag.) multiply-imaged system with three components, 
discovered during a surv ey for Extremely R ed Objects (EROs) in the fields of 10 massive galaxy 
cluster lenses at z ^ 0.2 (Smith et al.ll2002ar ). A redshift measurement on the brightest image, COa, 
was presented in lSmith et al.l ( 2002b ). Based on the 4000 A-break identification at A = 1.04 ± O.Ol/xm, 
the redshift is z = 1.60 ± 0.03. This image was also included in two LRIS masks, but we failed to 
detect strong spectral features in the wavelength range 5500-9200 A . 

Because of the clear symmetry of images COa and COb with respect to the critical line, it is possible to 
isolate 4 bright knots in each image a, b and c and use the corresponding 12 images to constrain the 
mass model. Furthermore, part of the south-west knot of COa is located within the radial caustic line 
in the source plane. We are able to detect a very faint radial arc predicted by the mass model (labeled 
CO-radial on Fig. [3]) after removing the light from the BCG. 

CI source 

The brightest cornponen t Cla of this system was observed during our NIRSPEC critical lines survey 
fsee lStark fc Ellis! I2OO6I ) and Stark et al. (2006b), in press). Three strong emission lines were identified 
as [O7JJ]5007, [0//7]4959 and H/3 (Figure g left), which unambiguously give a redshift z = 1.583. 
Optical spectroscopy of component Clc, included in one of our LRIS masks, could not identify any 
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strong [On] emission for this source, with a, 3 a upper hmit of 50 A for the equivalent width in rest- 
frame. 

• C4 source 

Although not multiply-imaged, this blue arc has very strong Lyman-a emission corresponding to z ~ 
2.63. An extended Lyman-a blob has also been detected in the IFU data (see Sect. 15.41 and Figure 
131). The mass model gives a very high magnification factor of /it ~ 4.0 mag., because the source is very 
close to a cusp. 

• C6 and C20 sources 

We identify a spatially-extended emission line in the LRIS spectrum at A = 4485 A , surrounding 
the C20c image (Fig. 21 middle). The lack of other strong emission lines in the optical range and 
the predicted position for two other components for C20 from the mass model confirms this to be 
Lyman-a at z — 2.689. In addition, this redshift is close to the Lyman-a blob associated with C4, and 
we identify a similar strong emission line at a similar redshift, z — 2.693, for the nearby source C25 
(~ 8" away in the source plane). 

The location and color of the lower surface brightness arc, C6, close to images C20a and C20b, is 
compatible with three merging images forming a single giant arc at a similar redshift (Fig. |4l middle). 
This is further evidence of a possible group of galaxies at z ~ 2.6 including sources C6, C20 and C25. 

• CI 5 source 

A strong asymmetric emission line is clearly seen on the LRIS spectrum at the position of CI 5a and 
C15b (FigurelH top-right), with a central peak at 7808 A. In both cases, a faint continuum is detected 
on the red side of the emission line, with a flux ~ 4.4 ± 1.5 x 10~^^ ergs s~^ cm~^ A~^. We interpret 
the emission line as Lyman-a at z — 5.421. This redshift could be slightly overestimated due to an 
unknown amount of self-absorption in the Lyman-a emission line and the absence of other spectral 
features. By averaging the extracted spectra of both images (Figure 01 top- right), the integrated flux 
measured within the emission line, of 9.7 ± 0.8 x lO"'^^ ergs s~^ cm~^, corresponds to a rest-frame 
equivalent width of ~ 35 A. We measure a Full- Width Half-Maximum of ~ 7 A when fitting the 
average emission line with a gaussian profile. The strong Lyman-a emission line of this source seems 
to dominate the overall /-band flux. 

• C23 source 

The identification of the pair of images C23a and C23b, predicted by the lensing model at z = 3.1, 
strengthens the Lyman-a interpretation for the single emission line seen in the spectrum of C23a at 
A = 5028 A. A fainter counterimage, C23c, predicted by the mass model, lies beyond the detection 
limit of the HST image. 

• C2 source 

This system is composed of two bright symmetrical images, identified close to the cluster center. Using 
the mass model, we predict a redshift of z ^ 1.4 for this source, with a magnification factor of fj, ^ 2.9 
mag. for C2a and C2b. A less magnified (/i ~ 1.8 mag.) counter-image C2c is also predicted by the 
model, but is not detected in the less-sensitive region of the HST image, at the junction between 2 WF 
chips. 

• CIO source 
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Fig. 4. — (Left) NIRSPEC bidimensional spectrum of component Cla, showing the detection of 3 strong 
features identified as [O///] and H/3 emission hnes at z — 1.583, both in positive (black) or negative (white), 
due to dithering offsets used. (Middle) Location of the LRIS slit and extended Lyman-a emission identified 
around image C20c (thick ellipse). Two brighter components C20a and C20b are predicted by the mass 
model, in agreement with the location of the critical line at z = 2.68 (solid curve). The shape and colors of 
the adjacent giant arc C6 are compatible with a similar redshift, suggesting a physical connection between 
these two sources. (Right) Upper panel: close-up of the 2-D sky-subtracted LRIS spectrum, showing both 
strong emission lines seen at the position of C15a and C15b. Lower panel: average extracted spectrum 
of C15, with the same wavelength range, revealing the clear asymmetrical shape of the spectral line. The 
relative sky background level is presented as a dotted curve. 
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A faint curved arc near a bright cluster member was identified by S05 as a multiple imaged candidate 
CIO. The optical colors of this blue arc are indeed compatible with a source at z < 2.5. It is probably 
formed by two images merging on the critical line; its shape is in agreement with predictions from the 
mass model. However, the low surface brightness of this arc and the lack of spectroscopic information 
strongly limit this interpretation, making it more uncertain than the systems previously mentioned. 
Therefore, we only include this source in Table [3] as a possible additional multiple imaged system. The 
mass model predicts a redshift of ~ 1.1 which is in very good agreement with the photometric redshift 
estimate of Z phnt = 1.15 =b 0.05 der ived from the broad band colors using the Hyperz photometric 
redshift code (jBolzonella et al.ll2000r ). A fainter counter-image (ClOc) is predicted by the mass model, 
but remains undetected in the HST image. 



Mass model 



4.1. Modeling method 



Using the new redshift measurements and the identification of further multiply- imaged systems, we are 
now in a position to considerably improve the precision of the mass model presented by S05. In doing so, 
we will maintain the Pseudo-Isothermal Elliptical Mass Distribution model (PIEMD, Kassiola & Kovner 
1993) adopted by S05 to infer the dark matter mass distrib ution. This parametric method has been used 
for modeling galaxy clusters as well as individual galaxies ( Covone et al. 2006 : Nataraian et al. IQQSh . It 
assumes each dark matter clump can be parameterized by a central position, ellipticity a/6, position angle 
0, central velocity dispersion ctq and two characteristic radii Vcore and rcut- The total mass of this profile is 
proportion al to rr„,t . A more d etailed discussion of the validity of this approach, in contrast to alternatives, 
is given bv lLimousin et al.l (|2007l ). 



The cluster galaxy population is incorporated into the len s model as galaxy-sca le perturbations to the 
cluster potential, assuming a scaling relatio n M/L = Const (IBrainerd et al.l ll996L S05) for all the Vcore, 
rcut and co- This is motivated by the similar Faber fc Jackson (Il976 ) scaling relation observed in elliptical 
galaxies. Following the same procedure as [Limousin et al. r(|2007l) . we keep the r*^^ and cTg values corre- 
sponding to a L* elliptical galaxy as free parameters, while keeping r*^^^ at 0.15 kpc. We select cluster 
galaxies in the field of Abell 68 by plotting the characteristic cluster red sequences {B-RJ02) and {Rjm-K) 
in two color-magnitude diagrams, keeping the objects pertaining to both red sequences. This reduces the 
photometric catalog used by S05, containing 69 galaxies with K band photometry, to 47 cluster members. 
This method is more efhcient than a single red sequence selection as we did not select any of the spec- 
troscopically con firmed background o r foreground ga laxies. As de scribed in S05, the K-hand photometry 
was obtained by Balogh et al. (2002) using GIM2D ( Simard 19981 ) to fit the surface brightness profiles of 
the cluster galaxies. This method gives more accurate results than one could obtain with SExtractor, be- 
cause SExtractor usually overestimate the local background around the brightest and most extended galaxies. 



In a first attempt at modeling the lens, we adopted the parameters given by S05, who included two dark 
matter halos. We found any attempt to reconstruct the mass distribution using a single clump was unable 
to reproduce the multiply- imaged systems accurately, confirming the strong bimodality of this cluster. By 
incorporating our spectroscopically-confirmed multiply-imaged systems, a reoptimization is necessary. To do 
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this, we use the new bayesian optimization method provided by the Lenstool softwaran (|Kneiblll993l Julio et 



al. in prep.) so that we can derive error estimates for each optimized parameter. The software optimizes the 
locations of each system in the source plane, based on the following estimator, which defines the goodness 
of the fit: 

i 

where x| is the same estimator for a given multiply- imaged system i, constructed by comparing the predicted 
positions of the N observed images in the source plane {xij,yi,j) to their barycenter {xf,yf): 

] = 1..N 

cFpos is the uncertainty in measuring the position of each image in the source plane. We use a typical value 
of (7/ = 0.2" for the same uncertainty in the image plane, and relate it to Upos with the amplification A: 

(7j = A CTpo^. 

For its important influence on the position of systems C15 and C20, we choose to optimize the velocity 
dispersion and the cut radius of the third brightest cluster galaxy. We also optimize the same parameters 
for the BCG and the galaxy #106 (adjacent to image COc), in order to match the location of all images from 
the CO system. 

The associated reduced is ^ 1-7, with 21 degrees of freedom, and the astrometric error on the position 
of the predicted multiple images is 0.37" in the image plane. The new parameters of the mass model are 
presented in Table ID The optimized values are slightly higher than the ones we obtain by using only their 
luminosity as a parameter for the gravitational potential. However, those values are not sufficiently high to 
propose a third dark matter clump. As additional confirmation for the quality of the mass model, we also 
checked that all spectroscopically-confirmed background sources for which wc did not identify multiple-image 
systems were predicted to be singly-imaged. 

In comparison with the results from S05, we have refined the mass model by increasing the number of 
constrained parameters from 11 to 21, while keeping a similar reduced value. In particular, we optimized 
the location of the second dark matter clump, where we have the main differences in the confirmed multiply 
imaged systems, and constrained individual parameters for two particular galaxies, in complement to the 
BCG. 



4.2. Results from the Mass Model 

Our improved mass model of Abell 68 enables us to compute a lensing mass of the cluster by integrating 
the derived surface mass density within a given projected distance R. Within a physical radius R < 500 
kpc, we obtain a value Mi^nsmg = 5.31 ± 0.17 x lO"'^'* M©. This is somewhat higher than the value of 
4.4±0.1 X 10"'^^ Mq derived by S05. By comparing both models using the individual parameters given in 
tableUl this difference arises mainly because of the higher mass of the second dark matter clump, for which we 
derive higher ctq rcut values compared to S05. The remaining model parameters (the location, ellipticity 
and orientation) lie within 3(7. We argue that the parameters of the secondary clump are better constrained 
in our model, given the identification of new multiple images to the NW of the cluster. The difference 



^For more information and to download the latest version of the code, see |http : //www. oamp. f r/cosmology/lenstool | 
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with S05 is further revealed by computing the ratio Mcen/Mtot between the mass of the central component 
(with sole contributions from the main clump and the BCG) and the total mass within 500 kpc. We find 
Mcen/Mtot — 0.56 ± 0.05, lower than the value of 0.68 derived by S05. This strengthens the bimodal nature 
of Abell 68. Our parameters for the individual contributing galaxies are quite similar to S05, except for the 
low r*„( value, indicating that the haloes have a smaller spatial extent. 

Within the uncertainties, we are now able to predict the location and expected fluxes of the counter- 
images for the observed multiple systems C2, CIO and C23. Furthermore, the model gives us an estimate of 
the redshift for the multiple systems C2 and CIO which do not yet have spectroscopic measurements. These 
predictions are summarized in Table [3] as bracketed values. 

In the process of the bayesian optimization, the software computes the magnification factors and the 
related error estimates, for each source included in the spectroscopic catalog or the multiple images catalog. 
These values are reported in Table [5] and [31 respectively. 

Apart from the multiple-image systems described before, all other background sources in our spectro- 
scopic catalog (C4, C7, C8, C12, C14, C24, C25, C26, C27) are predicted to be singly- imaged. This is 
compatible with our morphological data. For other sources that do not show multiple images (C3, C5, C9, 
Cll, C13, C18), we use the mass model to predict their maximum redshift z^ax- As the radius of the 
critical line increases with source redshift z^, a multiple image is expected if Zg > Zmax- These Zmax values 
are summarized in Table [3] and are consistent with the observed optical colors. Image C18 is predicted to 
be a single image at any redshift. 

5. Physical Properties of Faint Lyman-a Emitters 

Via our spectroscopic campaign and improved mass model, we are now in a position to explore the 
physical properties and luminosity distribution of a large sample of intrinsically faint 2 < z < 6 sources. We 
recognize that the volumes sampled may be unrepresentative but that our survey illustrates, by example, 
the promise of more extensive surveys that may soon be possible with larger cluster samples. 

First we selected from our spectroscopic catalog a subsample of seven high redshift sources [z > 1.5), 
which all have clearly detected Lyman-a emission lines. The measured rest-frame equivalent widths W 
exceed 2 A, and the majority have > 10 A . By contrast, typical LAEs selected within narrow-band 
surveys have > 30 A . The strong magnification thus gives us insight into the physical properties of these 
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[0.043] 


150±34 


179±7 


#106 


[-10.0] 


[-14.5] 


[1.2] 


[0.] 
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0.15 
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Table 4: Most likely parameters of the lens model with Icr error bars. From left to right: identification, 
astrometric position relative to the Brightest Cluster Galaxy, PIEMD parameters. 9 orientation is increasing 
from north through east. Bracketed values are not optimized. 



faint LAEs, including the star formation rate, the steUar mass and the physical scale corresponding to the 
star-forming regions. These physical values are derived and summarized in Table [51 



5.1. Star Formation Rates 



We explore two ways of determi ning the instanta neous star form ation rat e (SFR ) for a high redshift 
LAE, both using the calibration from Ecennicutt (1998) who assumed a IS aloetej \l9F>^ initial mass function 
with mass limits 0.1 and 100 Mq. The first calibration is based on the UV continuum luminosity iisoo: at 
1500 A in rest-frame, with the following relationship: 



5Fi?uv(MQyr-i) - 1.05 x lO^^o ^^^^^ (ergs s"! A 



(3) 



We estimated the individual iisoo values in this sample through the broad-band photometric measure- 
ments. 



The second calibration adopted from iKennicuttI ( 19981 ) is based on the intrinsic luminosity within the 
Lyman-a emission line, assuming no extinction and case B recombination (Brocklehurst 1971) : 



5FiZLya(M0yr-i) = 9.1 x lO^^^ i^y^lergss-i) 



(4) 



The corresponding SFR estimates are given in Table [51 

In most cases, we find good agreement between these two estimates, with a trend for 5'Fi?Ly,Q to be 
lower than SFRtjy, with a mean ratio SFR^y / SFRi^y a ~ 1 — 4. We interpret this difference as due to the 
specific properties of Lyman-a emission, which usually show some self-absorption or dust ext inction. The 
ratio above is quite sim ilar to that, ^ 3, typically found in high-redshift L yman- a sample s ( Santos et al 
2004; Aiiki et al. 2003 ). as well as for the most distant galaxies at z ~ 6.5 ( Hu et al. 2002t Kodaira et al 
20031 ). We notice two exceptions for C27 and C25; these show much higher SFR ratios (~ 50). Since C27 



shows a double-nucleus, its UV emission may be coming from an AGN. For C25, we estimate a 0.6"-offset 
between the center of the star-forming region and the center of the LRIS slit (originally aligned on C15a 
and C15b). Conceivably we missed the majority of the light coming from this object, which may explain the 
discrepancy. 



5.2. Stellar Masses 



Next we combine the constraints from multiband photometry and spectroscopic redshifts to derive 
estimates of the stellar mass associated with each source. A key question is whether LAEs are being seen at 
a special stage in their evolution, for example with a high star formation rate but low stellar mass. In this 
calculation, prior to SED fitting, it is important to remove the contribution of the Lyman-a flux from the 
photometric measurements. 

We derive the stellar mass using the Bayesian stellar mass code de veloped bvlBundv et al. which 
compares the SED of each object with a grid of synthetic SEDs from iBruzual fc CharlotI (|2003r ). assuming 
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a IChabried (j2003r ) Initial Mass Function. The star-formation history is parametrized as an exponential 
decaying burst SFR oc exp —t/r. In addition, the code gives an error estimate for each stellar mass, combining 
photometric errors and degeneracies in the model parameter space (age, reddening, metallicity, star-formation 
history). At 2 < z < 5.5, photometric measurements or upper limits in the 7 filters from Table I^TT] cover 
the restfi'ame UV and optical wavelengths, which is the main limiting factor in deriving stellar masses. 
Depending on the number of photometric bands (2 to 7) where a LAE is detected, the typical errors arising 
from the degeneracies in model parameters cover the range 0.1 to 0.5 in log]^Q(M*). It was not possible to 
derive a stellar mass for C26, because it is only detected in the Rf702W filter. 

Stellar masses, corrected for magnification, are reported in Table [H We also compute the invers e of the 



specific star formation rate (star formation rate per unit stellar mass. iBrinchmann et al.l (e.g. l2004J l) based 
on the SFR\jv estimate. This gives an indication on which timescale the star formation is taking place in 
each object. We find typical values of 500 Myr to 1 Gyr, suggesting that higher star formation may have 
occured in these LAEs during the past, unless they were formed very early. 



5.3. Physical scales 

In addition to brightening the observed flux of background sources, the magnification also stretches the 
angular sizes of the lensed images. This affects all the object shape parameters (a, b, 6), increasing the 
observed solid angle by the same factor fj,. We are thus able to probe a physical scale R in the source plane, 
with: 

B=^Dj,{z) (5) 

where Da{z) is the angular distance for a source at z. 

For most of the high redshift sources which are not resolved (along one direction) in the HST image, this 
value of R is an upper limit (Table O. In most cases, we find that we are able to probe high redshift sources 
at sub-kiloparsec scales. Another interesting physical property we derive is the intrinsic surface density of 
star formation (S = SFRuy / (ir i?^)), which is a lower limit when sources are unresolved. 



5.4. Lyman-a emission surrounding C4 

We now turn to a specific issue relating to the origin of the so-called Lyman-a blobs, of which C4 at 
z — 2.63 is the first lensed example. We have used the updated mass model to reconstruct the source-plane 
morphology of C4 . The high magnification (^ 4 mags.) enables us to resolve the source morphology at ~ 30 
pc scales. The galaxy displays a bright component with an intrinsic elongated shape and several knots, at 
a scale length of ^3 kpc along the major axis. A fainter component having a similar shape and one bright 
knot is located 4.6 kpc away (Fig. [51 fourth panel). Such an elongated shape and small physical scale 
is not uncommon among high redshift galaxies: Ravindranath et al. ( 20061 ) have measured that 50% of 
galaxies in their sample of z > 2.5 LBGs taken from GOODS show a bar-like morphology, and scale lengths 
of about 1.7 — 2. kpc. 



The strong Lynian-a emission of C4 detected with LRIS was also identified in the VIMOS/IFU data 
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C27 


C4 


C20c 


C25 


C23a 


C26 


C15ab 


Ellis ot al. 


z 


1.75 


2.63 


2.69 


2.69 


3.13 


3.68 


5.42 


5.58 


fi (mags) 


1.72 


4.15 


3.61 


1.12 


2.37 


2.06 


2.79 


3.80 


iisoo^ 


6.1 


5.5 


0.41 


4.6 


0.39 


0.24 


3.68 


0.21 




4.6±0.6 


9.3±0.6 


4.8±0.4 


1.3±0.5 


1.5±0.3 


2.6±0.2 


9.7±0.8 


6.8±0.7 




2.1±0.3 


1.6±0.1 


1.4±0.1 


2.9±1.1 


1.4±0.28 


3.4±0.3 


29±2.4 


6.9±0.7 




6.4 


5.8 


0.44 


4.8 


0.41 


0.25 


3.9 


0.22 




0.19±0.03 


0.15±0.01 


0.13±0.01 


0.26±0.09 


0.13±0.03 


0.31±0.03 


2.62±0.22 


0.35±0.04 


logu-)iM*/MQ) 


9.6 ±0.1 


8.7 ±0.2 


8.4 ±0.3 


9.7 ±0.1 


8.6 ±0.5 


N/A 


9.4 ±0.3 


8.2 ±0.2 


iJ(kpc) 


2.8 


2.2 


< 0.36 


< 0.72 


< 0.37 


< 0.34 


< 0.28 


< 0.19 




0.25 


0.38 


> 1.08 


> 2.94 


> 0.95 


> 0.68 


> 15.8 


> 1.93 




6.2 


0.9 


5.7 


10.4 


9.7 


N/A 


6.4 


7.2 



Table 5: Summary of physical properties derived for all LAEs at z > 1.5 in our sample, and comparison 
with the mean values of the source a,t z ^ 5.6 found bv lEllis et al.l (j2001l ) in Abell 2218 (rightmost column). 
From top to bottom: redshift, magnification factor fi, unlensed UV luminosity {L\) at rest-frame A =1500 
A, observed integrated flux in the Lyman-a emission line, corresponding unlensed Lyman-a luminosity, star 
formation rate derived from the UV continuum (SFRuy) or the Lyman-a emission (SFRi^y^), stellar mass 
estimate from SED fitting (see Sect. 15. 2p . intrinsic physical scale R (see Sect. 15. 3p . star formation surface 
density S and specific star formation rate. 



h 17 1 — 2 

X 10 ergs s cm 
'^xlO'*-'^ ergss"-*^ 

f X 10* yr 



.1 





^ Inwffflplpjii] 






\ 






■IT 



Fig. 5. — (Left) Negative color image of a 40"side region around the source C4, produced by combining three 
different wavelength slices from the IFU datacube : a 15 A-wide region encompassing the Lyman-a emission 
(blue), and two broad wavelength regions at redder wavelengths (green and red). The extended emission 
around the arc C4 is clearly detected in the blue frame. 

(Right) Lyman-a emission properties of C4 observed with the IFU data. From left to right: rest-frame 
equivalent-width mapping of the Lyman-a emission (with corresponding color-scale) , contours of Lyman-a 
emission overplotted on the HST image, i?-band reconstruction of C4 in the source plane. 
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(Fig. [5l first panel). Compared with the observed size of the arc in the continuum image R702 (4x1.2") 
this emission is significantly more extended (~15x3", third panel of Fig. [5]). After correction for an average 
magnification, this corresponds to an intrinsic scale of 10 kpc. 

We use the i?702W image to estimate the continuum fiux of C4 at the wavelength of the Lyman- 
a emission, after applying a fc-correction factor assuming a UV spectral slope of 2.0 typical of starburst 
galaxies. By comparing the spatial coverage of the Lyman-a emission and the continuum, we construct 
a map of rest-frame Lyman-a equivalent width Wq (Fig. [51 second and third panel). The Wq values are 
typically in the range 5 < Wq < 10 A in the region detected on the HST image. However, the map shows 
a significant increase of Lyman-a equivalent width, with values Wq > 10 A, in the northern part of the 
arc, where the continuum is hardly visible in R band. Because of the size of the blob and the presence 
of this stronger stellar continuum at the center of C4, the most probable mechanism producing such an 
extended Lyman-a em i ssion is the pr e sence of a superwind outflow originating from the central starburst 



Taniguchi et al.l (|2001f) : IWilman et all (|2005f ) 



We also observed C4 with 2 different LRIS configurations, located across or along the longer dimension 
of the arc (Fig. [5]). In both cases, we observed two different components to the Lyman-a line on the 2D 
spectrum : a stronger emission centered around 4411 A and a fainter emission region at a slightly bluer 
wavelength (4408 A ). The LRIS data also shows a tilt in the emission region, due to an offset in the 
central peak of the Lyman-a line (top panel of Fig. [6]). These components may be related to the different 
star-forming regions identified morphologically in the HST image. Unfortunately, because of the poorer 
resolution, no similar offset in the central wavelength or variations in the line profile were detected in the 
IFU data. 

The main difference between C4 and the very luminous Lyman-a blobs is its very sm all physical size 



2.2 kpc, compared with the i? > 30 kpc selection criteria for Lyman-a blobs adopted bv iMatsuda et al 



(2004), along with the rest-frame equivalent width criterion Wo > 20 A . For this reason, it is more likely 
that the associated physical processes may be different, the size of C4 being more similar to extended Lyman- 
a emssions observed around star-forming galaxies. In this strongly lensed case, the effects of outflows on 
the Lyman-a emission can be studied in higher details by way of the lensing magnification, stretching the 
observed scales. 



5.5. Intrinsic properties of the z ^ 5.4 source 



We finally turn to the most distant lensed sour ce, C15, at z = 5. 4, and whether it is a further example 



of the intriguing lensed source detected at z—5.6 bv lEUis et al.l (j200ir ) in AbcU 2218. 

C15 has an intrinsic Lyman-a luminosity of ^ 3 x 10^^ ergs s~^. T his is about 5 times fainter than 
typic al values found for z ~ 5.7 LAEs targetted by narrow-band searches (IHu et al.ll2004i IShirna-saku et al. 
2006 ). but not too dissimilar to those in the much fainter z ~ 6.5 sample from Kashikawa et al.l "( 2006 ) who 
used very deep (~ 10 ks) spectroscopy with Keck and Subaru to confirm 17 faint LAEs. 

All three images of this source are unresolved along the shear direction, and therefore suggest a very 
small physical scale R < 300 pc in the source plane. More interestingly, we infer a minimum star formation 
surface density of ^ 16MQkpc~^, which is much higher than for any other object in our sample, making 
this source similar to the most active star forming regions at any redshift. 



The strong Lyman-a emission dominates the optical fiux and in this respect it is very similar to the z ~ 
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5.6 object found by I Ellis et al.l (|20011 ) in the cluster Abell 2218. To compare these sources, we have performed 
a re-analysis of the photom etric data for that object, taking into account HST data that arrived for Abell 
2218 subsequent to the EUis et all J2OO1I) an alysi s. This includes HS T/ACS, NICMOS and Spitzer/IRAC 
observations presented by 



Kneibetal 



|2004h and lEgami et all |2005[ ). We find that the EUis et al source is 
now detected in the Ifsuw, zpssoLP, Jpiiow and ffieoiv filters (Fig. [7] and Table [6]). The source is faint 
(AB~ 26.6) and remains undetected with IRAC at 3.6 and 4.5 microns. However, the IRAC upper limits 
are quite high (around AB~24), due to the proximity of the BCG (^ 15") and coarse point spread function. 

In order to compare the revised data with that for C15, we compute and incorporate in the last column of 
Table [5] various physical properties for the Ellis et al. object alongside those for the other Lyman-a emitters. 
When deriving the stellar mass, we note the significant change implied by the new data (M* ~ 10^ M© c.f. 
1O^~^M0 given by Ellis et al). The major change arises from different assumptions : in order to reproduce 
the observed Lyman-a emission and the non-dete ction of the UV spect ral continuum in LRIS, Ellis et al 
used the STARBURST99 spectrophotometric code lLeitherer et al.l (jl999f ) to derive an upper limit of 2 Myrs 
on the age of the source, assuming a constant SFR and no extinction. 

Although the new HST photometry of this object is restricted to the rest-frame wavelength A < 2500 
A, the SED fitting method used in our reanalysis should be more reliable because it is based on fewer 
assumptions, especially on the star-formation history. Our best SED fitting model also predicts Spitzer/IRAC 
continuum fluxes consistent with the non-detection of this pair in the 3.6 iim and 4.5 fim bands. 



Compared with the z ~ 5.6 galaxy from lEllis et al.l (j200ll ). C15 is less magnified and is intrinsically 
more luminous and more massive, although with a similar physical size. The mass model predicts the source 
to be located ~ 2.4 kpc from the caustic line in the source plane. This suggests we are actually seeing a 
true small isolated object, since any similarly bright region close to C15 would also be highly magnified and 
multiply imaged. 

Another difference between the two sources arises from the rest-frame UV stellar continuum, which 
is detected in the LRIS spectrum of the z ^ 5.4 source redwar d of the Lyman-g emission. Running the 
STARBURST99 code on this object in a similar way as done bv^ Ellis et al. (2001), the upper limit on the 
age is much larger, typically 50 — 100 Myrs, with a constant SFR of 2.6 Mo.yr"-'^. This gives a total mass of 
1 — 3 X 10* Mq, which is closer to the SED fitting estimate, in comparison with the source in Abell 2218. 



At the time of its discovery, the lEUis et al.l (j200ll ) object revealed a very small and low-mass (10^~^) 
source at z ~ 5.6, which was also reported to be very young. With the new photometric reanalysis, and 
by finding another example of a LAE with similar size and Lyman-a flux, it is more probable that we are 
observing two lO'* to 10^ Mq objects with very different star-formation histories. The Ellis et al. source 
previously formed the majority of its stellar mass and is having a very young burst, whereas C15 has been 
forming stars for a larger time scale. 



Vfqo6W If814W 



ZF850LP 



Jf 



now 



Hf160W 



IRACa.e^im IRAC4.5^m 



> 28.14 26.94 ±0.18 26.67 ±0.15 26.60 ± 0.28 26.62 ± 0.33 > 24.1 



> 24.3 



Table 6: Broad-band photometry of the z = 5.6 galaxy discovered bv lEUis et al.l (|200lh . measured on the 
new HST-ACS/NICMOS and Spitzer/IRAC images. 
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Fig. 6. — LRIS observations of C4. Left panels : closeups of each 2D spectrum showing the two components 
appearing in the Lyman-a emission. Right panels : corresponding LRIS slit configuration plotted over the 
HST image. 




Fig. 7.— HST broad-band detections of the z ^ 5.6 pair found bv IeIHs et all |200lh in Abell 2218. From 
left to right: /f814H'j ^ssolp, ^/fiiow, Hfi^wv- North is up and East is left. The light from the nearby 
BCG (located towards the lower-right corner) has been subtracted from the NICMOS data for clarity. 
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6. Discussion 



We have demonstrated in our survey of Abell 68 that, after applying magnification correction, we can 
reach LAEs with typical unlensed absolute AB magnitudes ~ —20 < A/ab < —16.5 in the rest-frame UV. 
Th is is about 2 magnit udes fainter than the UV selection criterion for typical Lyman Break Galaxies used 
by ISteidel et al.l (l2003r) . of R\n < 25.5 at z ^ 3 and the faintest LAEs found with narrow-band searches 
(e.g., Fvnbo et al. 2003t Gawiser et al. 20061 : Shimasaku et al. 20061). Moreover, the lens ed Lyman-a sources 
have much lower equivalent widths, similar to those seen in LBGs (|Shaplev et al.ll2003l ). 



The derived properties of the Lyman-a emitters are dependent on the individual magnification factors 
fj, used to correct each image for the lensing effects. The error estimates on fj,, computed by the bayesian 
optimization method, range from 1 to 16 % and were quadratically added to the measurement errors in 
Table [5l in the case of Li^y a , SFRi^y „ and Af * . In each case, the photometry is the dominant source of 
error. We acknowledge these values of /d are still dependent on the approach used to parametrize the mass 
model, for instance the use of the PIEMD profile for the dark matter haloes. Nevertheless, the Lyman-a 
emitters presented here are strongly lensed and have been identified in the central regions of the cluster, 
where the magnification factors are associated with the location of the critical lines. These critical lines are 
constrained by the same set of multiple images, quite independently of the parametrization used for the mass 
model (PIEMD or Navarro- Frenk- White (NEW) profile). Therefore, we are confident that the computed 
magnification factors are reliable. 

Although such a lensing survey is not strictly flux-limited, we can use this unique probe to gauge the 
properties of the faintest LAEs yet located at z > 2. We find the typical stellar masses are log in(Af<,/Af(Ti) ^ 
8.5 to 9.5. These values are similar to the typical stellar masses of the bright LAEs from iGawiser et al 



(j2006l ). found by selecting very high equivalent width {Wq > 150) LAEs at z=3.1. Our LAEs have fainter 
UV and Lyman-a luminosities, and therefore a lower SFR, by typically 1 to 2 magnitudes. Even if the 
wavelength range covered by the broad-band photometry dominates the errors in the stellar masses derived 
in both surveys, this is indicative that, at comparable stellar masses, our LAEs are more quiescent than the 
objects found in narrow-band searches. 

Recognizing the limitations of our sample, as a point of illustration, we compare on Fig. [8]the luminosity 
range of our lensed emitters with current constraints on the cumulative Lyman-a luminosity function at 
z ^ 3, mostly based on narrow-band searches. The number of dedicated searches for LAEs at z 3 from the 
literature is quite limited. We report in this diagram the compilation of 5 samples, either in blank fields or 
overdense environments, presented by Fvnbo et al. (2 0?). as well as an estimate of the number density from 
a Lya survey in the GOODS-south field (jNilsson et al. (2006) and Nilsson et al. 2007, to be submitted). 
Based on our subsample of 6 emitters at 2 < z < 3.7 and the surface area probed in the high magnification 
region (/j, > 0.75 mag.) of the source plane, we estimate a number density of 4.31^ 7 ^ 10~^ Mpc~^ down to 
L = 1.4 X 10*^ ergs s~^. Because our spectroscopy only sparsely covered this region, we report this value as 
a lower limit in Fig. [H Although we acknowledge the limits due to very small statistics and the small areas 
probed in these lensed surveys, our sample of z '--^ 3 emitters gives constraints at fainter luminosities that 
any of the narrow-band searches at this redshift. More observations are however needed in order to better 
constrain the faint-end slope of this luminosity function. 

An interesting question is the li kelihood of finding highly magnified sources at z ^ 5.5 like C15 and the 
source reported by I Ellis et al.l (|200l[ l : we can estimate this through the space density of i?-dropouts in our 
magnified field of view. Based on our photometric RIz color-color selection technique (see Sect. 13. ip . we 
identified one source in the field of view covered by the HST image. Since it is detected at ^ 8 cr level in the 
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Fig. 8. — Cumulat ive luminosity function of Lyman-a emitters. Red datapoints are z ~ 3 estimates from 



Fvnbo et al.l (|200ir ). either in overdense regions (open circles) or blank fields (crosses). We also give recent 
results at z ~ 3.1 from a Lya survey in GOODS-Sou th (Nilsson et al. 20 07, to be submitted). Black points 



Santos et al 



( 20041 ) and black curves are more recent 



correspond to the compilation of 2 ~ 5 surv eys from 
fits to the z ^ 5.7 Luminosity Function by Shimasaku et al.l (|2006l ). We overplot in blue number density 



estimates based on the sample of 6 Lya emitters at 1.7 < z < 3.7. 
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/-band filter, where i?-dropouts are brighter, and has a magnification factor of fi ^ 2.5, we computed the 
effective covolume in the source plane having /i > 0.75, taking into account the reduction of the surface area 
due to lensing effects. We derive an effective covolume of 420 Mpc"^, in the redshift range 5 < z < 6 probed 
by our set of filters. Assuming Poisson noise statistics, we obtain a rough estimate of 2.4l|^ xlQ-^Mpc"^ 
for the number counts of i?-dropouts. Although cosmic variance effects can be quite large when probing such 
a small volume, this result is in good agreement with the space density of low-luminosity Lyman-a sources 
at z ~ 5, on the above-mentioned diagram (Fig. [8]). 

From these simple calculations, C15 does not seem to be a serendipitous case of low-luminosity z ^ 5 
LAE, since we would expect to find one such source in the magnified region of the cluster. By observing 
a larger number of lensed fields, such as the new HST clusters imaged by ACS, we could build a more 
significant sample of similar objects and compare their physical properties. 



7. Conclusions 

We have performed a spectroscopic analysis of background sources in the field of the massive cluster 
Abell 68 and identified 26 lensed images in the range 0.3 < z < 5.5, including 7 Lyman-a emitters at z > 1.7. 
Using the new redshift measurements and identification of new multiple-image systems, we perform a precise 
modeling of the cluster mass distribution with 5 spectroscopic systems, and also predict the redshift and 
counter-images for two remaining systems. This makes Abell 68 one of the best-modelled lensing clusters, 
along with Abell 1689 and Abell 2218, allowing for the precise measurement of its dark- matter distribution. 

We derived the star formation rates, stellar masses and physical scales for our sample of high-redshift 
Lyman-a emitters. The broad-band luminosities of these objects is comparable to the faint LAEs found in 
deep narrow-band searches, but their equivalent widths are much lower, making them 1 to 2 magnitudes 
fainter in Lyman-a luminosity. Two of these sources show a more extended Lyman-a emission region, when 
compared to the stellar continuum in the HST image. For one of them, we demonstrated how we use IFU 
data to probe regions with distinct Lyman-a equivalent widths. The stretch provided by lensing enables us 
to characterize small regions of Lyman-a emission which would otherwise be beyond reach of ground based 
Integral Field Instruments. Although the large equivalent widths are comparable to giant Lyman-a blobs 
observed around massive forming galaxies, we interpret the small physical scales of these lensed emissions 
as outflows originating from a central starburst. 

The highest redshift z 5.4 multiple-im age source of this sample is very similar to the pair of strongly- 



lensed images identified bv I Ellis et al.l (|200lh in Abell 2218, in terms of magnification and physical size, albeit 
being intrinsically more massive and more luminous. W e therefore expect to detect it in IRAC images of 
similar depth as the data presented in lEgami et al. ( 2005[ l. in the first two channels of this instrument. Such 



measurements would tighten the constraints on the stellar mass of this source, by reducing the degeneracies 
in the other model parameters, and also provide additional information on its age and star-formation history, 
when combined with the properties of the Lyman-a emission. 

Our survey provides the first tentative indications of the density of faint Lyman-a emitters at z ^ 3, 
down to unlensed fluxes of ~ 2 x 10^^. Although our survey is not fiux- limited nor complete in any formal 
sense, the cumulative Lyman-a luminosity function we derive indicates the promise of a dedicated search 
for lensed Lyman-a emitters through larger samples of well-mapped clusters now being surveyed with HST, 
and the ability of this approach to complement narrow-band searches carried out in blank fields. 
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